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, Charge transport through a short DNA oligomer (Dickerson dodecamer) in presence of structural 
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fluctuations is investigated using a hybrid computational methodology based on a combination of 
quantum mechanical electronic structure calculations and classical molecular dynamics simulations 
with a model Hamiltonian approach. Based on a fragment orbital description, the DNA electronic 



structure can be coarse-grained in a very efficient way. The influence of dynamical fluctuations 
arising either from the solvent fluctuations or from base-pair vibrational modes can be taken into 
account in a straightforward way through time series of the effective DNA electronic parameters, 
evaluated at snapshots along the MD trajectory. We show that charge transport can be promoted 
through the coupling to solvent fluctuations, which gate the onsite energies along the DNA wire. 
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I. INTRODUCTION 



The electrical response of DNA oligomers to applied voltages is a highly topical issue which 
has attracted the attention of scientists belongig to different research communities. The 
variability of experimental results is reflected in DNA being predicted to be an insulator, 1 
a semiconductor, 2 ' 3 or a metallic-like system. 4 ' 5 This fact hints not only at the difficulties 
encountered in carrying out well-controlled single-molecule experiments, but also at the dra- 
matic sensitivity of charge migration through DNA molecules to intrinsic, system-related or 
extrinsic, set-up mediated factors: the specific base-pair sequence, internal vibrational ex- 
citations, solvent fluctuations, and the electrode-molecule interface topology, among others. 
As a result, the theoretical modelling of DNA quantum transport remains a very challeng- 
ing issue that has been approached from many different sides, see e.g., Refs. 6,7,8,9 for 
recent reviews. While most of the models originally used started from a static picture of 
the DNA structure, 10 ' 11 ' 12 ' 13 it has become meanwhile clearer that charge migration through 
DNA oligomers attached to electrodes may only be understood in the context of a dy- 
namical approach. 11,15,16,17,18,19,20,21 Hole transfer experiments 22 ' 23 ' 24 ' 25 ' 26 had already hinted 
at the strong influence of DNA structural fluctuations in supporting or hindering charge 
propagation. Hence, it seems natural to expect that dynamical effects would also play a 
determining role in charge transport processes for molecules contacted by electrodes. A 
realistic inclusion of the influence of dynamical effects onto the transport properties can 
however only be achieved via hybrid methodologies combining a reliable description of the 
bio molecular dynamics and electronic structure with quantum transport calculations. Ah 
initio calculations for static biomoecular structures can provide a very valuable starting 
point for the parametrization of model Hamiltonians; 27 ' 28 ' 29,30 ' 31 ' 32 ' 33 ' 34 ' 35 however, a full first- 
principle treatment of both dynamics and electronic structure lies outside the capabilities 
of state-of-the-art methodologies. 

In this paper, we will ellaborate on a recent study on homogeneous DNA sequences 2 " by 
addressing in detail some methodological issues. Our focus will be on the so called Dicker- 
son dodecamer 36 with the sequence 3 — GCGCTTAACGGC— 5' and for which the effect of 
the dynamical fluctuations becomes very clear. Our approach combines classical molecular 
dynamics (MD) simulations with electronic structure calculations to provide a realistic start- 
ing point for the description of the influence of structural fluctuations onto the electronic 
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structure of a biomolecule. Further, the information drawn from such calculations be used 
to formulate low- dimensional effective Hamiltonians to describe charge transport. A central 
point is the use of the concept of fragment orbitals 37 ' 38 ' 39 ' 40 which allows for a very efficient 
and flexible mapping of the electronic structure of a complex system onto a much simpler 
effective model. Taking a single base pair as a fragment and considering only one fragment 
orbital per base pair, we end up in a tight-binding Hamiltonian for a linear chain where both 
onsite energies €j(t) and electronic coupling terms Vjj+i(t) are time-dependent variables: 



The dynamical information provided in this way builds the starting point of our treatment 
of quantum transport through biomolecular wires. In the next Section, we briefly describe 
the computational methodology used to obtain the effective electronic parameters of the 
model Hamiltonian, which will be then introduced in Sec. Ill, where we also illustrate how 
to relate the Hamiltonian of Eq. (1) to a different model describing the coupling of a time- 
independent electronic system to a bosonic bath. Further, expressions for the electrical 
current as well as the relation between the auto-correlation function of the onsite energy 
fluctuations and the spectral density of the bosonic bath will be derived. Finally, we discuss 
in Sec. IV the transport properties of the Dickerson dodecamer in vacuum and in presence 
of a solvent. We stress that in contrast to other models which explicitly contain the coupling 
to vibrational excitations or to an environment 41 ' 42 ' 43 ' 44 ' 45 at the price of introducing several 
free parameters, our methodology potentially contains the full dynamical complexity of the 
biomolecule as obtained from the MD simulations. One main advantage of our approach 
is the possibilty to progressively improve the degree of coarse-graining by an appropriate 
re-definition of the molecular fragments. 

II. MOLECULAR DYNAMICS AND MODEL HAMILTONIAN FORMULATION 
A. Computational methodology 

We will first give an overview of the fragment-orbital approach used in our computations; 
further details can be found elsewhere. 38 ' 39 ' 46 

The core of our method is based on a combination of a charge self-consistent density- 
functional parametrized tight-binding approach (SCC-DFTB) and the fragment orbital 




(1) 
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FIG. 1: Left panel: Schematic representation of the fragment orbital method used to perform a 
coarse-graining of the DNA electronic structure. A fragment consists of a single base pair (not 
including the sugar phosphate backbones). As explained in the text, the hopping matrix elements 
Vj j+i between nearest-neighbor fragments are compued using the molecular orbital basis of the 
isolated base pairs. These calculations are then carried out at snapshots along the molecular 
dynamics trajectory hence leading to time dependent electronic structure parameters. By keeping 
only one relevant orbital per fragment, the electronic structure can be mapped onto that of a 
linear chain (right panel). Transport observables can be computed at each simulation time step. 
Alternatively, the time dependence o the electronic structure (related to structural and solvent 
fluctuations) can be transferred to a bosonic bath as done in this paper. 

concept ; 46 both are used to compute the electronic parameters ej(t) and Vjj + i(t) for the 
effective tight binding model introduced in Eq. (1) in a very efficient way, see Fig. 1 for a 
schematic representation. The Vjj+x are calculated using the highest occupied molecular 
orbital $; computed for isolated bases as Vjj+i = (<&j\H\<&j + i) , where the $j's can be ex- 
panded in a valence atomic orbital basis rj^ on a given fragment: $j = Xlup^V The 
are obtained from calculations on isolated bases and stored for subsequent use to calculate 
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Vjj+\. Therefore, one can write: 

w = EE c X +1 ^i^>- ( 2 ) 

fj, V 

The Hamilton matrix in the atomic orbital basis H^ u = (ri^H^J) evaluated using the SCC- 
DFTB Hamiltonian matrix is pre-calculated and stored, thus making this step extremely 
efficient, i.e., it can be calculated for geometry snapshots generated by a classical molecular 
dynamics simulation even for several nanoseconds. Additionally, the minimal LCAO basis set 
used in the standard SCC-DFTB code has been optimized for the calculation of the hopping 
matrix elements, 38 ' 39 and the results are in very good agreement with other approaches. 28 ' 46 
Concerning now the coupling to the solvent, a hybrid quantum mechanics/molecular me- 
chanics (QM/MM) approach has been used, implemented in the SCC-DFTB code, 1 ' leading 
to the following Hamiltonian matrix in the valence atomic orbital basis rj v : 




S x " u A 

Here, Qs are the Mulliken charges of the quantum-mechanical region and the Qa are the 
charges in the MM region (backbones, counter-ions, and water), S^ v is the atomic orbital 



overlap matrix, H^ C ~ DFTB is the corresponding zero-order Hamiltonian matrix, and R a $ is 



the distance between the DNA atom where the AO orbital rj^ is located and the MM atom 
in the solvent with charge Qs- The last term explicitly takes into account the coupling to 
the environment (solvent). The €j = Vjj and Vjj±i from Eq. (2) using Eq. (3) can now 
be calculated along the MD trajectories 21 ' 39 . The off diagonal matrix elements strongly 
depend on structural fluctuations of the DNA base pairs, but they are only weakly affected 
by the solvent dynamics, the opposite holding for the onsite energies which are considerably 
modified by solvent fluctuations. 21 ' 39 The Fourier transform of the onsite energies auto- 
correlation function provides information about the spectral ranges which are more strongly 
contributing to the fluctuations of the electronic parameters, see also the next sections. 
Thus, we have found that the apparently most important modes are located around 1600 
cm 1 corresponding to a base skeleton mode, and at 800 cm -1 related to the water modes. 39 
Both contributions modulate the onsite energies significantly on a short time scale and a 
long time scale of about 1 ps, respectively. 
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B. Model Hamiltonian for electronic transport 



Using Eq. (1) directly for quantum transport calculations may mask to some degree different 
contributions (solvent, base dynamics) to charge propagation through the DNA 7r-stack. 
Moreover, since Eq. (f ) contains random variables through the time series, we are confronted 
with the problem of dealing with charge transport in an stochastic Hamiltonian. This is a 
more complex task which has been addressed e.g., in the context of exciton transport 48 ' 49 ' 50 ' 51 
but also to some degree in electron transfer theories. 52 ' 53 ' 54 Here, we adopt a different point 
of view and formulate a model Hamiltonian, where the relevant electronic system, in this 
case the fragment orbital-derived effective DNA electronic system, is coupled explicitly to a 
bosonic bath. The latter will encode through its spectral density the dynamical information 
drawn from the MD simulations on internal base dynamics as well as solvent fluctuations. 
The Hamiltonian can be written in the following way: 

H = H (, ', (l > - E ( ^ + ''•'••) ( 4 ) 

3 j 

—bath ~t~ -^tunnel -^lcads 

-f^bath = E ^-aB^Ba 

a 

#el-bath = E ^gd]dj(B a + B^) 

oi,j 

^ tunnel = ^ (* ks J C L d J + nX -) 
H\ eSL ds = ^ ^ C^sC^^Cks 

The time averages (over the corresponding time series) of the electronic parameters (ej) t and 
{Vjj + i) t have been split off to provide a zero-order Hamiltonian which contains dynamical 
effects in a mean-field-like level. The effect of the fluctuations around these averages is 
hidden in the vibrational bath, which is assumed to be a collection of a large (N — > oo) 
number of harmonic oscillators in thermal equilibrium at temperature k^T. The bath will 
be characterized by a spectral density J(w) which can also be extracted from the MD 
simulations, as shown below and in Ref. 55 . Since we are interested in calculating the electrical 
response of the system, the charge-bath model has to also include the coupling of the system 
to electronic reservoirs (electrodes). The coupling to the electrodes will be treated in a 
standard way, using a tunneling Hamiltonian H tunnc \ which describes the coupling to the s- 
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lead with s— left (L) or right (R). Later on, the so called wide-band limit will be introduced 
(the corresponding electrode self-energies are purely imaginary and energy- independent), 
thus reducing the electrode-DNA coupling to a single parameter. 

The previous model relies on some basic assumptions that can be substantiated by the 
results of the MD simulations: 39 (i) The complex DNA dynamics can be well mimic within 
the harmonic approximation by using a continuous vibrational spectrum; (ii) The simulations 
show that the local onsite energy fluctuations are much stronger in presence of a solvent than 
those of the electronic hopping integrals (see also the end of the previous section), so that 
we assume that the bath is coupled only diagonally to the charge density fluctuations; (hi) 
Fluctuations on different sites display rather similar statistical properties, so that the charge- 
bath coupling A Q is taken to be independent of the site j. This latter approximation can be 
lifted by introducing additional site- nonlocal spectral densities Jjj + i(u); this however would 
make the theory more involved and less transparent. In Fig. 2 we show typical normalized 
auto-correlation functions of the onsite energy fluctuations for the Dickerson dodecamer in 
both solvent and vacuum conditions as well as one case of off-diagonal correlations between 
nearest-neighbor site energies (inset). Also shown are fits to stretched exponentials, which 
are in general equivalent to a sum of simple exponential functions. This suggests the presence 
of different time scales and the non-trivial time-dependence of the fluctuation dynamics. 
The fits become obviously less accurate at long times due to the reduced number of sampled 
data points with increasing time (the oscillatory behavior becomes stronger). From the 
figure we first see that the off-diagonal correlations decay on a much shorter time scale as 
the local ones, so that on a first approximation their neglection can be justified; further the 
decay of the correlations for the vacuum simulations is considerably much faster than in 
a solvent indicating the strong influence of the latter in gating the electronic structure of 
the biomolecule. The corresponding correlation functions for the hopping integrals Vjj + i(t) 
(not shown) display even shorter relaxation times, so that the approximation Vjj+i(t) = 
(Vj,j+i(t)) t is enough for our purposes (the hopping integrals are self-averaging). 
In order to deal with the previous model, we first perform a polaron transformation of the 
Hamiltonian Eq. (4), using the generator U = exp[^^ a g a d\d(,{B a — B a )], which is nothing 
else as a shift operator of the harmonic oscillators equilibrium positions. The parameter 
9a = ^a/^a gives an effective measure of the electron- vibron coupling strength. Since we 
will work in the wide-band limit for the electrode self-energies, the renormalization of the 
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FIG. 2: Normalized auto-correlation function 5 ' 3 Ci(t) = (5ei(t)Sei(0)) / (Sef) and averaged nearest- 
neighbor correlation function C^i+ift) = (5ei(t)5ei + i(0)) / (5ej5ej + i) (inset) of the onsite energy- 
fluctuations. The solid lines are fits to stretched exponentials which suggests the existence of 
different time scales, a typical situation in the dynamics of bio-molecules. On average, the decay 
of Cj i i+i(t) occurs on a much shorter time scale than that of Ci(t), so that our model will only 
include in a first approximation local fluctuations, see Eq. (1). 

tunneling Hamiltonian by a vibronic operator will be neglected. As a result, we obtain a 
Hamiltonian with decoupled electronic and vibronic parts and where the onsite energies are 
shifted as (ej) t — > (ej) t — / °° doj J {uj) / uj . However, as it is well-known, 55 the retarded Green 
function of the system is now an entangled electronic-vibronic object that can not be treated 
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exactly; we thus decouple it in the approximate way: 



.41,57 



g nm (t,t') = -\e{t-t')^[d n {t)x\t),dlit')x{t')\ + ) (5) 

« -i9(t - t') {(d n (t)dl(f)) (x\t)X(t>)) + (dl(t')d n (t)) {X{t')X\t))} 
= 9{t - \G> m (t,t')e~^ - G< m (t,t')e-^' 



= J2^ 2 + N ^ e ~ iUat + N - e+iUat ] 

a 

In this equation, 0(£ — t') is the Heaviside function and the pure bosonic operator X(t) = 
ex PE a 9a(Ba ~ E a )]. in the last row of Eq. (5) we can pass to the continuum limit and 
express <p(t) in terms of the bath spectral density J(u;): 58 

. > i , , fvjj ,„ , . i r 00 , j(w) . 

(pit) = — / au; — — coth- (1 — cosu;t) — l — / du — -sinwt. (6) 

n Jo u k B T hj u 2 

C. The electrical current 



We derive in this section the expression we are going to use to calculate the electrical current 
through the DNA oligomer under study. Starting point is the well-known Meir-Wingreen 
expression for the current from lead s: 59 

Is = 2 ij ^ Ti ^(e)g>(e)-j:>(e)g<(e)}- 

Now, we can exploit the decoupling approximation used in Eq. (5) together with the wide- 
band limit in the electrode-molecule coupling to write e.g., for the left electrode: 

X<(E)G>(E) = j ^f L (E)(l-f R (E>))t(E>ME-E>), 

/j -pi 
^(l-h(E))ME')t(E'ME'-E), 

$(£) = J ^ e ^ e -*W 

Hereby we have used the explicit expressions for the greater- and lesser-Green functions: 

G>'<(E) = G (E)(^<(E) + £>'<(£))Gto, (7) 

the index indicating that the vibrational degrees of freedom have already been decou- 
pled. The transmission-like function t{E) is given by t(E) = Tr |Go(£ , )rLGj(i?)rR |. The 
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retarded matrix Green function Gq(E) is calculated without electron-bath coupling but in- 
cluding the interaction with the electrodes: Gq (E) = E + irj — H + iT L + iT R , H being 
the electronic part of the Hamiltonian of Eq. (4). Using these results, the right-going current 
can be written as 



Tj^J ^ t{E ' ] {f^m-fR(E')ME-E>)-(l-f L (E))f R (E>)$(E'-E)} 



a similar expression holding for the left-going current, when the indices L and R are inter- 
changed. The total current can be written in a symmetrized way: Jt = (Jl — -^r)/2. By 
looking at Eq. (6), two limiting cases can immediately be obtained: the zero charge-bath 
coupling (<f) = 0) which implies <&(E) = 2tt5(E). In this limit we recover the conventional 
expression for coherent transport, involving only the transmission function t(E). In the 
high-temperature and/or strong coupling limit to the bath, a short-time expansion of <fi(t) 
can be performed, yielding: 

t 2 f°° , T/ x , fa ,t f, J(v) 

W = 7TT / dujJ{uj)coth——-i- / duj , 

2hJ k B T hj uj 

E 

"'therm'' 1 <" 

n 

In the former expression ^^thcrm is related to an inverse decoherence time. In the high 
temperature limit K t herm ~ k B TE reorg . The Fourier transform of the previous expression can 
be calculated straightforward and gives: 



$(E) = a — exp 



7T . (E + E, 



\2 



reorg) 



crm 

Thus, the current calculated using Eq. (8) becomes a convolution of t(E) with a Gaussian 
function. Notice the similarity of this expression with a Frank-Condon factor appearing in 
the Marcus electron transfer theory. 



D. Spectral density from molecular dynamics 

The next issue at stake is how to relate the bath spectral density J(uj) to the time series 
of the electronic parameters as obtained through the molecular dynamic simulations. To 
illustrate this point we will consider the simple case of a single level, whose site energy is a 
Gaussian random variable, coupled to left and right electrodes according to the Hamiltonian 
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(we assume for simplicity that (e(t)) = 0): 



H = 5e(t)d)d + H tmmcl + H lcads (9) 

To deal with this problem we may use equation of motion techniques for the retarded dot 
Green function G(t,t') = — (i /K)6(t — t') ({d(t), c^(t')}) in the time domain: 

(ihdt - Se{t))G(t,t') = S(t - t') +J2 t la G Ut,t'), (10) 

with Gk a (t,t') = —{i/h)6{t — t') ({c kjQ (t), d^(t')}Y A similar equation of motion for the 
latter function allows to introduce the electrode self-energies: 

(i hd t - 5e(t))G(t, t') = 5(t -t') + J drE(t, t)G{t, t'), 

£(t,T) = 5> kia | 2 e- ie -* (11) 

Note that the Green function thus obtained is still a random function, since no average over 
the random variable 5e(t) has been yet performed. Using the wide-band limit in the lead 
self-energies £ a (i, £') = — iT a 5(t — t') we get, with r = Y L + Y R : 

(i hd t - Se(t) + i r)G(t, t') = 5{t - t'). (12) 

Introducing the time-exponential U(t,t') = exp(— (i/ft) f*,ds(Se(t) — iT) allows to write a 
closed solution for the dot's Green function: 

G(t,t') = - 1 -9(t-t')U(t,t') (13) 

The average Green function over the random variable Se(t) yields now in energy-space: 

• poo roc 

(G{E)) = -- dte L ^ t - t 'HU(t-t')) = -i / dte^ E+ir X t -*Ue-yt> dsS < s A. (14) 

60 



Performing now a cumulant expansion'' of the averaged exponential and taking into account 
that cumulants higher than the second one exactly vanish due to the Gaussian nature of the 
fluctuations and to the fact that Se(t) is a classical variable, we get (with if — 0): 



POO 

( G (E)) = - 1 - I dte^ E+ ^ 1 ti ds ^ ds ' ^ {s)&e{s ' ] \ (15) 



o 



which is the formally exact solution of the problem. We may now look at the same problem 
from a different point of view by considering explicitly the coupling of a single site with 
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time-independent onsite energy to a continuum of vibrational excitations. Using the polaron 
transformation together with the approximations introduced at the beginning of Sec. II B, 
we can write the retarded Green function as: 

G(E) = -- / dteM E+ir *e-*®, (16) 
Jo 

where <f>(t) has been already defined in Eq. (6). By comparison of Eqs.( 15) and (16), there 
must exist a relation between the (real) correlation function and the real part of 4>(t). The 
latter can be re-written in the following way: 

Reom — — duo — —coth- (1 — cosa;t) 

H J uo 2 k B T 

— ds ds' < — I duo J(uo) coth - — — cos [uo (s — s')] > , 
Jo Jo I " Jo ^b-T J 

so that from this we conclude that 

(5e(s)5e(s / )) = h [ duo J(uo) coth -p^- cos [uo (s — s')l. 
Jo k sT 

From here it follows by inversion: 

J(o;) = — tanh- — — / dt cosuot C{t) = — tanh - — —j(uo). (17) 
Tin k^T J 7rn k^T 

The previous result was obtained in Ref. 55 in the weak coupling, perturbative regime to the 
vibrational system; we have shown here that it can be also extended to the case of arbitrary 
coupling. Similar expressions could be derived for (spatially) non-local spectral densities. 



E. Results 

To illustrate our methodology we have focused on the Dickerson dodecamer (DD) which has 
a non- homogeneous base sequence. In contrast to our previous study 20 where homogeneous 
DNA sequences were addressed like poly(G)-poly(C) and poly(A)-poly(T), the study of the 
DD nicely illustrates the role of the solvent fluctuations in gating the electronic structure. To 
study the effect of the conformational dynamics as well as the fluctuations of environment 
onto the charge transport properties, we have performed classical MD simulations using 
the AMBER-parm99 force field92 with the parmBSCO extension93 as implemented in the 
GROMACS 1 software package. The static geometries were built with the 3DNA program 62 
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FIG. 3: Upper panel: Time averaged transmission function for simulations performed in solvent 
and in vacuum. Lower panel: Coherence parameter C{E) for the Dickerson dodecamer in solvent 
and in vacuum. Notice that the influence of the solvent fluctuations is to spread out the spectral 
support of C{E). As expected, the coherence parameter is on average larger for the vacuum 
case where a large part of the fluctuations is suppressed due to the absence of an environment. 
Nevertheless the transmission in the vacuum case is much smaller, nicely illustrating the positive 
influence of the environment in gating charge migration. 
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FIG. 4: Time averaged of the absolute value of the onsite energies along the segment covering 
bases 3 and 10 of the Dickerson dodecamer, i.e., over the sequence 3 -GCGCTTAACGGC-5 ) 
for simulations in vacuum (upper panel) and in solvent (lower panel). The outer most two bases on 
each end were not included in the calculations to avoid undesired boundary effects. The averaged 
energy profile in presence of the solvent becomes smoother but also the fluctuations around the 
averages are stronger. This smoothing reduces the energy barriers between the sites and hence 
favours charge migration. 

while the starting structures for the MD simulations were created using the make-na server. 63 
After a standard heating procedure followed by a 1 ns equilibration phase, we performed 30 
ns MD simulations with a time step of 2 fs. The simulations were carried out in a rectangular 
box with periodic boundary conditions and filled with 5500 TIP3P water molecules and 22 
sodium counterions for neutralization. Snapshots of the molecular structures were saved 
every 1 ps, for which the charge transfer parameters were calculated with the methodology 
described in the previous sections. When speaking of simulations in vacuum conditions, 
we mean that the last term in Eq. 3 is left out. Upon mapping the electronic structure 
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Voltage [V] 

FIG. 5: Electrical current for seven base pairs of poly(dG)-poly(dC) and poly(dA)-poly(dT) 
oligomers in vacuum. The current is considerably suppressed when comparing with the case where 
solvent fluctuations are included. 

onto a linear chain as discussed in the former sections, we have first computed the time 
and energy dependent quantum-mechanical transmission function T(E,t) by evaluating for 
every set of charge transfer parameters, i.e., at each simulation time step t, (T(E,t)) t = 
(l/TyiD)^2jT(E,tj). This quantity is expected to provide some qualitative insight into 
different factors affecting transport (solvent vs. vacuum) but its use is restricted by the 
fact that for longer chains inelastic, fluctuation mediated channels will play an increasing 
role and then a pure elastic transmission can not catch all the transport physics. In this 
latter case, the model Eq. (5) seems to us more appropriate since it includes the dressing 
of the electronic degrees of freedom by the structural fluctuations. In Fig. 3 we show the 
averaged transmission function (upper panel) for the cases where solvent effects are included 
or neglected, respectively. We clearly see that solvent-mediated fluctuations considerably 
increase and broaden the transmission spectrum. Its fragmented structure is simply due to 
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the fluctuations in the onsite energies, which make the system effectively highly disordered. 
Another way of looking at the effect of fluctuations is via the introduction of a coherence 
parameter 20 which is defined as C(E) = (T(E, t)) t / (T(E, t) 2 ) t . If the fluctuations are weak 
C(E goes to one, while a strong fluctuating system will lead to a considerably reduction 
of C(E). Obviously, this parameter has only a clear meaning over the spectral support 
of the transmission function. C(E) is shown in the lower panel of Fig. 3; the coherence 
parameter in vacuum conditions can become larger than in solvent for some small energy 
regions due to the reduction of the fluctuations. However, this does not turn to be enough 
to increase the transport efficiency for the special base sequence of the Dickerson DNA, 
since this system has in the static limit a distribution of energetic barriers due to the base 
pair sequence. To further illustrate the influence of the solvent, we have plotted in Fig. 4 
the time averaged onsite energies along the model tight-binding chain. Remarkably, the 
presence of the solvent " smoothes" the averaged energy profile (though the amplitude of the 
fluctuations clearly becomes stronger). To compute the current through the system, we use 
the Hamiltonian Eq. (5) and the current expression, Eq. (8), derived with this model. The 
Fermi energy in these calculations was fixed at the upper edge of the transmission spectrum, 
see Fig. 3. The qualitative results are however not considerably changed by slightly changing 
its position. In Fig. 5 the current is shown for the two cases of interest. Due to the presence 
of tunnel barriers in the wire which are not fully compensated on average by the gating 
effect of the environment, the absolute current values are rather small when compared with 
those of homogeneous sequences-". However, the current including the solvent is roughly 
fifteen times larger than for that obtained from the simulations in vacuum. Though our 
model Hamiltonian in Eq. (5) does not fully contain all the dynamical correlations encoded 
in the time dependent electronic parameters, we nevertheless expect that their inclusion 
would lead to an even further increase of the difference between solvent and vacuum results, 
thus suppporting our main conclusions. 



III. CONCLUSIONS 



We have presented a hybrid methodology which allows for a very efficient coarse-graining of 
the electronic structure of complex biomolecular systems as well as to include the influence of 
structural fluctuations into the electronic parameters. The time series obtained in this way 
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allow to map the problem onto an effective low dimensional model Hamiltonian describing 
the interaction of charges with a bosonic bath, which comprises the dynamical fluctuations. 
The possibility to parametrize the bath spectral density J(oS) using the information obtained 
from the time series makes our approach very efficient, since we do not need to use the 
typical phenomenological Ansatze (ohmic, Lorentzian, etc) to describe the bath dynamics. 
The example presented here, the Dickerson dodecamer, shows in a very clear way that 
solvent-mediated gating may be a very efficient mechanism in supporting charge transport 
if the static DNA reference structure already posseses a disordered energy profile (due to 
the base sequence). Our method is obviously not limited to the treatment of DNA but it 
can equally well be applied to deal with charge migration in other complex systems like 
molecular organic crystals or polymers, where charge dynamics and coupling to fluctuating 
environments plays an important ro le. 52 ' 53 ' 54 ' 64 ' 65 ' 66 ' 67 ' 68 
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